Flavor Physics and CP Violation Conference, Vancouver, 2006 



Hadronic B Decays in PQCD 

Satoshi Mishima 

School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, U.S.A. 



I review recent progress on exclusive hadronic B meson decays in the pcrturbative QCD approach, with focus on 
puzzles in the branching ratios and the CP asymmetries of the B — > ttK and B —* tttt modes, and polarization 
fractions in B — > VV modes. 



1. Introduction 

B factory experiments have accumulated a lot of 
data and have reported many interesting results [1]. 
Some observables, mixing-induced CP asymmetries 
for b — > s penguin modes, branching ratios and di- 
rect CP asymmetries for B — > ttK and tttt, and po- 
larization fractions for penguin-dominated B — > VV 
modes, have exhibited some deviations from naive ex- 
pectations in the Standard Model. 

It is necessary to go beyond naive estimations for 
understanding the observed deviations. The pertur- 
bative QCD (PQCD) approach [2, 3] is one of the 
theoretical attempts to include subdominant contri- 
butions, such as spectator and annihilation diagrams, 
and higher-order corrections. PQCD has applied to 
various two-body B decays at leading order (LO) in a s 
and has made reasonable predictions for various decay 
modes. Recently, the important next-to-leading-order 
(NLO) contributions were evaluated in the B — > ttK, 
tttt, and pp decays to investigate the discrepancies be- 
tween the LO PQCD predictions and the data [4, 5]. 

This talk is organized as follows: In Sec. 2, I briefly 
review the PQCD factorization formula. I discuss 
the branching ratios and the CP asymmetries of the 
B — > ttK and 7T7T decays with the NLO corrections 
in Sec. 3. The LO PQCD predictions of the polariza- 
tion fractions are presented in Sec. 4. Section 5 is a 
summary. 



2. PQCD Factorization Theorem 

Most of the calculations of B decay amplitudes rely 
on the factorization of decay amplitudes into a prod- 
uct of short-distance and long-distance physics. QCD- 
improved factorization (QCDF) [6] and soft-collinear 
effective theory (SCET) [7] are based on collinear fac- 
torization theorem, but PQCD is based on kr factor- 
ization theorem. 

Employing collinear factorization theorem, some 
decay amplitudes involve a singularity arising from 
the end-point region of parton momentum fractions. 
An end-point singularity implies that a decay ampli- 
tude is dominated by soft dynamics and cannot be fac- 
torized. Such soft contributions are regarded as phe- 
nomcnological parameters, which can be fitted from 
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Figure 1: Factorization of hadronic two-body B meson 
decays in the PQCD approach. 

the experimental data. 

In the PQCD approach with kr factorization the- 
orem, the Sudakov factor ensures the absence of the 
end-point singularities [8] . All amplitudes can be fac- 
torizable into parton distribution amplitudes <J>, the 
Sudakov factors e~ s , and a hard kernel H: 

A(B -> M 2 M 3 ) = $m 2 ® $m 3 ® H <g> e" s ® $ s , (1) 

where ® stands for convolutions in both longitudi- 
nal and transverse momenta of partons [9]. The 
schematic picture of the factorization theorem is given 
in Fig. 1. The distribution amplitudes, which are uni- 
versal in the processes under consideration, are de- 
termined from experiments, the light-cone QCD sum 
rules, lattice calculations, or other theoretical meth- 
ods. The hard kernel is characterized by a hard scale 
Q ~ y/Krn,b, where A is a hadronic scale and rrib the b 
quark mass [10-12], and can be evaluated as an expan- 
sion in powers of a s (Q) and A/Q. The hard kernels of 
the spectator and annihilation contributions, as well 
as the emission contribution, are calculable and start 
from 0(a s ). PQCD predicts a large direct CP asym- 
metry in B° — > tt^K^ as a result of a large strong 
phase arising from annihilation penguin diagrams [2] . 



3. B — > ttK and Puzzles 

The current data of the direct CP asymmetries of 
B —> ttK and the branching ratios of B — > 7T7t [1], 

A C p(B ± -> tt ^) = (4 ±4)%, 



fpcp06_423 



2 



Flavor Physics and CP Violation Conference, Vancouver, 2006 



A CP (B° - 
Br(B° - 
Br(P° 



n^K*) = (-10.8 ±1.7)% , 
> ti^tt^) = (4.9 ±0.4) x 10~ 6 , 
-►ttV) = (1.45 ± 0.29) x 10~ 6 , (2) 



have been inconsistent with the expected relations 

A CP (B ± - tt ^) « A CP (B° -> 7r^X±) , 
Br(P° -» t^tt*) > Br(P° -» 7r°7r°) . (3) 

These relations can be understood in the topological- 
amplitude decompositions [13]: 



A(B+ - 


- 7T+if°) 


V2A(B+ - 


- 7T if+) 


A(B° - 


> 7T-if+) 


V2A(B° - 


-> 7r°X°) 



p' 



_pl _ pi 



ew 



P' 



P' 



(T' + C) e^ 3 



(4) 



and 



A{B° - 
V2A(B+ - 
V2^(P° 



7T 7T J 

► ir+ir°) 
- ttV) 



-T — Pe*^ 2 
-T 1 C P ew e" 
-C+(P-P ew )e^ , (5) 



o«<P2 



where TW, CW, PW, and P e ( ^ stand for the color- 
allowed tree, color-suppressed tree, penguin, and 
electroweak penguin amplitudes, respectively, and 
<f>2 and <f>3 are the weak phase defined by V u b = 
\Vub\exp(-i<f>3), V td = \V t d\ exp(-i^i), and <j> 2 = 
180° — <f>i — <t>3- Assuming the hierarchies, P' > 
T", P^ > C" for P -» ttK and T > C, P > P eu) 
for P — > 7T7T [14-16], the relations in Eq. (3) can 
be derived. The current data seem to require a 
large electroweak penguin amplitude [16-22], a large 
color-suppressed tree amplitude [23-26], or both of 
them [27-29]. 

The LO PQCD predictions follow the naive expec- 
tations as shown in Tables I and II [2, 3, 30]. Because 



Table I Branching ratios for the B 
in units of 10 -6 . 



nK and nir decays 



Mode 


Data [1] 


LO [30] 


NLO [4, 5] 


B ± - 

B ± - t^K* 

B° V^K* 

B° - tt ^ 


24.1 ± 1.3 
12.1 ±0.8 
18.9 ±0.7 
11.5 ± 1.0 


14.4 - 26.3 
7.9 ~ 14.2 
12.7 ~ 19.3 
4.5 ~ 8.1 


23.6l 1 8 4 i 
13.61 1 ";? 
20.4i 1 ^ 

8.7^ 1° 


_B° -» t^tt* 
S± - 77*77° 
B° -> 77°77° 


4.9 ±0.4 
5.5 ±0.6 
1.45 ±0.29 


5.9 ~ 11.0 
2.7 - 4.8 
0.10 - 0.65 


6-51 U 

4.ol tt 
o-29±8:S 



LO PQCD predicts a negligible C { '\ NLO contribu- 
tions may enhance sufficiently. Recently, the most 
important NLO contributions from the vertex correc- 
tions, the quark loops, and the magnetic penguin were 
calculated in the PQCD approach [4, 5]. The sum of 



Table II Direct CP asymmetries for the B 
decays in percentage. 



nK and nn 



Mode 


Data [1] 


LO [30] 


NLO [4, 5] 


B ± -> n ± K° 


-2 ±4 


-1.5 0.6 


0±0 


B ± - tt ^ ± 


4±4 


-17.3 10.0 




B° -> yr^A-i 


-10.8 ±1.7 


-21.9 - -12.9 




B° -> n°K° 


2 ±13 


-1.03 0.90 


7+3 
' -4 


B° - 77*77* 


37 ± 10 


16.0 ~ 30.0 


1 8+ 20 


S± - 77*77° 


1 ±6 


0.0 


0±0 


5° - 77°77° 


9 o+40 
^°-39 


20.0 - 40.0 


63±|S 



the latter two reduces the LO penguin amplitudes by 
about 10% in the B — > wK decays and affects the CP 
asymmetries little. The vertex corrections enhance 
C by a factor of three and rotate the phase of C , so 
that V ± C is parallel to P' ± P' ew . Thus, the direct 
CP asymmetry of B ± — » 7r°i4T ± vanishes as shown in 
Table II. 

To compare the predicted B — > ttK branching ratios 
with the data, it is useful to consider the ratios 



R 

R c 

Rn 



Br(P° -^tt*K ± )t b + 
Br(P± ir ± K°) t b o 
Br(P± -» 77° if*) 
Br(P± -» 7r±if°) 

1 Br(P° -» ti^X*) 

2 Br(P° 77°^°) 



= 0.85 ±0.06 , 

1.00 ±0.08 , 

0.82 ±0.08. (6) 



As shown in Fig. 2, the PQCD predictions for R and 
P c are consistent with the data, whereas that for R n 
is larger than the data [4] . The larger R n is due to the 
smaller P° — ► tt K branching ratio, which might im- 
ply a significant P' ew . However, the discrepancy is not 
serious at this moment, because the theoretical un- 
certainty is representative and the actual uncertainty 
could be larger. 

The enhanced C also affects the mixing-induced 
CP asymmetry of B° — ► tt°Ks, given by 



S w 



2Im(A 7r o Ks ) 

1+ lA^Ksl 2 



with 



K°K S = ~ e 



pi _ p. 



it 

ew 



P' 



P' 

ew 



(7) 



(8) 



S^oks may deviate from S cSs — sin(2</>i) in b — > ccs 
modes. Therefore, it is essential to include the NLO 
corrections in the calculation of S n ox s ■ It was found 
that the predicted possible deviation of S 7r a Ks is still 
small, S^ Ks -S c5s = 0.06±g;g| [4]. Moreover, PQCD 
predicts a positive deviation, opposite to the observed 
one. 
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Figure 2: R, R c , and R n as functions of c/>3 from NLO PQCD with the bands representing the theoretical 
uncertainty [4]. The two dashed lines denote la bounds from the data. 



C for B — > 7T7T, as well as C for _B — > ttK, is 
enhanced by the vertex corrections, but it is insuf- 
ficient to accommodate the B° — > tt°tt° branching 
ratio to the measured value [4]. NLO PQCD pre- 
dicts \C/T\ w 0.2 for B -> tttt, though a much 
larger |C/T| w 0.8 is required to explain the observed 
B° — > 7r 7T branching ratio [25]. 

The same NLO PQCD formalism was applied to the 
B — > ,0/9 decays [5], which are sensitive to the color- 
suppressed tree contribution. The predicted B — > pp 
branching ratios are listed in Table III. The LO results 



Table III Branching ratios for the B — > pp decays in 
units of 10" 6 [5]. 



Mode 


BABAR [1] 


Belle [1] 


LO 


NLO 


B° -» pV 
B ± ^P ± P° 
B° - p°p° 


30 ± 4 ± 5 
17.2 ±2.5 ±2.8 
< 1.1 


22.8 ± 3.8+J;;? 
31.7±7.lt|;? 


27.8 
13.7 
0.33 


nc q+25.3 
13.8 

i6.o +1 8 5 ;? 

Q9+ 110 



differ from those in the previous LO analyses [31, 32] 
slightly due to the different choices of the hard scale 
and parameters. The NLO PQCD predictions for the 
B° — > p^p and B — > P^P° branching ratios are 
consistent with the data. Because the decay ampli- 
tudes for B — > pp are similar to those for B — > irn, 
the branching ratio of B° — > p°p° is expected to be 
larger than that of B° — > 7r°7r° due to the meson decay 
constants f p >f n - In fact, the NLO predictions follow 
this expectation, and the central value of the predicted 
B° — > branching ratio has almost reached the 

experimental upper bound. The NLO PQCD analysis 
has thus confirmed that it is unlikely to accommodate 
both the B° — > 7r°7r° and B° — > p°p branching ratios 
to the data simultaneously. Hence, the B — > 7T7t puz- 
zle is confirmed in the PQCD approach. All proposed 
resolutions to the B — > tttt puzzle should survive the 
constraints from the B° — ► data. 

I comment onB^ ttK and 7T7t results in other the- 
oretical approaches. In QCDF, is enhanced by the 
NLO jet function obtained from SCET and the large 
B° — » tt°tt° branching ratio can be explained [34]. 
However, the inclusion of the NLO jet function over- 



shoots the B° — > p°p° branching ratio and dete- 
riorates the predictions for the B ± — > tt^K^ and 
B° -► Tr^Xi direct CP asymmetries [5]. In SCET, 
incalculable soft contributions are regarded as phe- 
nomcnological parameters, which can be fitted from 
the experimental data. It was found that the charm- 
ing penguin, which is one of the phenomenological pa- 
rameters, is large as laW/rWl « 1 [35-37]. Conse- 
quently, the large B — > 7r°7r° branching ratio can be 
realized. The B° — ► p°p° branching ratio should be 
checked in the same formalism. The ratio jT^' is, 
however, real in the leading-power SCET formalism, 
and therefore the B ± -> tt°K ± and B° -> n^K* di- 
rect CP asymmetries can not be explained at the same 
time [35-37]. 



4. Polarizations in B — > VV 

In the naive factorization approximation, the longi- 
tudinal and transverse polarization fractions, Rl and 
i?y _l, respectively, in B — > VV modes obey the power- 
counting rules [38] , 

R L ~ l-0{m v /m 2 B ) , R\ U ~ 0(m 2 v /m 2 B ) , (9) 

where mg is the mass of the B meson and my that 
of the emitted vector meson from the weak vertex. 
The data of the longitudinal polarization fractions are 
given by [1] 



Rl(B° 




— u.au 1 _o.028 ' 


Rl(B+ 


-pV) 


= 0.96 ±0.06 , 


Rl(B + - 


► K*+p°) 


- 91+ 23 

— U.J1_ 21 , 


Rl(B + - 


> K*°p+) 


_ n 4s+ 09 

- U - 48 -0.08 > 


Rl(B+ - 


-> <f)K* + ) 


= 0.50 ±0.07, 


Rl(B° 




= 0.48 ±0.04. 



The polarization fractions for tree-dominated modes, 
B — > pp, satisfy Eq. (9). However, those for penguin- 
dominated modes, B — > cf>K* and K*°p + , obviously 
conflict with the naive expectations. 
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For the penguin-dominated modes, the polarization 
fractions could be modified by sub-leading contribu- 
tions. The penguin annihilation contribution from the 
(S - P)(S + P) operators, which follows R L ~ R\\ ~ 
R±, could decrease the longitudinal fraction [38]. In 
LO PQCD, the spectator and penguin annihilation 
contributions help to reduce Rl as shown in Table IV, 
but it is not enough to explain the B — ► <f>K* and 
K*°p+ data [39-41]. R L for B+ -> K*+p° remains 



Table IV Polarization fractions in B — > cjiA'*, (I) without 
spectator and annihilation contributions, and (II) with 
spectator and annihilation contributions [39]. 



Mode 






«ll 


R± 


0|| (rad) 


0_l (rad) 


B° -> <f>K*° 


(I) 


0.923 


0.040 


0.035 


7T 


7T 




(II) 


0.750 


0.135 


0.115 


2.55 


2.54 


B+ -> ^A"* 4 - 


(I) 


0.923 


0.040 


0.035 


7T 


7T 




(II) 


0.748 


0.133 


0.111 


2.55 


2.54 



as Rl ~ 0.85, which is consistent with the data, be- 
cause this process involves additional tree amplitudes. 
The tree-dominated modes, which are insensitive to 
the sub-leading corrections, follow the naive counting 
rules in Eq. (9) [5, 31]. Therefore, only the penguin- 
dominated modes B — ► <pK* and K*°p + have exhib- 
ited anomalies in the measured polarization fractions. 

There have been several mechanisms proposed to 
explain the observed B — > cf)K* polarizations. In 
PQCD, it was proposed that the B — > K* form fac- 
tor Ao, associated with the longitudinal polarization, 
may be smaller than the central value of the LO 
PQCD prediction [42]. Postulating a smaller value, 
Ao w 0.3, which does not contradict to any exist- 
ing measurements, Rl for B — ► (f>K* decreases to 
0.6. Another mechanism in PQCD was proposed in 
Ref. [32]. Adopting a modified definition of the hard 
scale, which is a source of theoretical uncertainty, Rl 
for B — > <j)K* could approach to 0.6. However, Rl for 
B + — > K*°p + becomes about 0.8, which is inconsis- 
tent with the data. This is because the sign of the real 
part of the annihilation amplitude for B + — > K*°p + 
is opposite to that for B — ► [32]. Small ifo for 

the penguin-dominated modes might come from the 
complicated QCD dynamics, but it is important to 
explain both the B — ► </>if * and data. 



5. Summary 

In this talk, I have summarized the recent works 
on exclusive hadronic B meson decays in the PQCD 
approach, concentrating on the observed deviations in 
the branching ratios and the CP asymmetries of the 
B — > wK, 7T7T modes, and the polarization fractions of 
penguin-dominated B — ► VV modes, which are sensi- 
tive to subdominant contributions. 



Including the important NLO contributions, the 
color-suppressed tree amplitude is enhanced by the 
vertex corrections, and therefore the predicted direct 
CP asymmetries of the B — ► wK modes become con- 
sistent with the experimental data. However, it is 
unlikely to accommodate both the B° — > 7r°7r° and 
B° — ► p°p° branching ratios to the measured ones si- 
multaneously. 

The polarization fractions of the penguin- 
dominated B — ► VV modes deviate from the 
naive power-counting rules, including the spectator 
and annihilation contributions. However, it is not 
enough to explain the observed data. A small 
longitudinal fraction for B — > (f>K* might come from 
QCD uncertainty, but it is necessary to explain both 
the B — » (j)K* and K*°p + data. NLO corrections to 
the polarization fractions should be studied in future 
work. 
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